Magnetic field barriers in graphene: an analytically solvable model 



Enrique MilpatB and Manuel TorrcsQ 

Instituto de Fisica, Universidad Nacional Autonoma de Mexico, 
Apartado Postal 20-364, Mexico Distrito Federal 01000, Mexico 

Gabriela Murgui'a@ 
Facultad de Ciencias, Universidad Nacional Autonoma de Mexico, 
Apartado Postal 21-092, Mexico Distrito Federal 04021, Mexico 
(Dated: January 13, 2013) 

We study the dynamics of carriers in graphene subjected to an inhomogeneous magnetic field. For 
a magnetic field with an hyperbolic profile the corresponding Dirac equation can be analyzed within 
the formalism of supersymmetric quantum mechanics, and leads to an exactly solvable model. We 
study in detail the bound spectra. For a narrow barrier the spectra is characterized by a few bands, 
except for the zero energy level that remains degenerated. As, the width of the barrier increases 
we can track the bands evolution into the degenerated Landau levels. In the scattering regime a 
simple analytical formula is obtained for the transmission coefficient, this result allow us to identify 
the resonant conditions at which the barrier becomes transparent. 
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I. INTRODUCTION 

The discovery of graphene^—, a single layer of carbon in a honeycomb lattice has generated a lot of excitement, 
due to its unique electronic properties^— and its potential application in electronic devices. Electrons in graphene 
are described by a massless two dimensional relativistic Dirac equation 7 - - — , that yields a gaplcss linear spectrum near 
to the K and K' points of the first Brillouin zone. Graphene exhibits a variety of pseudo relativistic phenomena, 
providing an unexpected connection between condensed matter physics and quantum-relativistic phenomena. Among 
others we can cite: the Zitterbewegung and its relation with the minimal electrical conductivity at vanishing carrier 
concentration^*^, the unconventional quantum hall effect^iii, and the Klein tunneling^—. The Klein tunneling 
effect has important implications for the future design of graphene based electronic devices, because massless Dirac 
fermions cannot be effectively confined by electrostatic barriers, in particular for normal incidence the barrier becomes 
complete transparent^. 

Some schemes have been proposed in order to avoid the obstacle that represents the Klein tunneling, in order to 
confine electrons in graphene based structures. An interest proposal refers to the use of inhomogeneous magnetic fields 
that produce magnetic barriers^. Previous studies have consider the cases of single^, double^ or multiple magnetic 
barrier o 17 ' 18 . All of these refer to square magnetic barriers with sharp edges. In this paper we consider a magnetic 
barrier in which the edges are smoothed out. We select a magnetic field with an hyperbolic profile. We show that 
the corresponding Dirac equations can be analyzed within the formalism of supersymmetric quantum mechanics, and 
leads to an exactly solvable model. We study in detail the bound spectra. For a narrow barrier the spectra displays a 
series of bands separated by gaps, as the width of the barrier increases the bands evolve into the degenerated Landau 
levels. In the scattering regime a simple analytical formula is obtained for the transmission coefficient, this result 
allow us to identify the resonant conditions at which the barrier becomes transparent. 

The use of inhomogeneous magnetic fields has received considerable attention both in the experimental 2 ^— and 
theoretica l 20 ' 25 " — study of two-dimensional electron gases (2DEG) in semiconductor heterostructures. Various config- 
urations of local inhomogeneous magnetic fields have been created and studied, using microfabricated ferromagnetic 
and superconducting structures deposited on top of a 2DEG. Interesting transports phenomena have observed, among 
others: magnctoresistance and commensurable oscillations, anomalous transport along special e.g., snakelikc trajecto- 
ries, etc. Although there exist to date no experimental realization of similar configurations in graphene, they should 
be produced in the near future. The configuration for the hyperbolic magnetic field, see below Eq. ([1]), provides a 
good approximation to the shape of the magnetic barrier produced by a ferromagnetic film deposited in a 2DEG^i. 
Hence we expect, that apart of its intrinsic theoretical interest, the results obtained in this work will be useful in 
order to analyze the confinement by magnetic barriers in graphene samples. 

The paper is organized as follows. In section II we study the model for Dirac fermion dynamics in graphene when 
the system is subjected to an inhomogeneous magnetic field. We show that the system can be analyzed within the 
formalism of supersymmetric quantum mechanics. The effective potentials and the explicit analytical solution for 
the wave function are discussed. In section III we study in detail the bound spectra and analyze it corresponding 
degeneracy, both for a narrow barrier, and also in the limit in which the barrier width becomes comparable to the size 



2 



of the system. In section IV the dispersion regime is analyzed in detail. In section V the conclusions are presented. 



II. GRAPHENE IN AN INHOMOGENEOUS MAGNETIC FIELD 



We focus on an electron in a single graphene layer subject to an inhomogeneous perpendicular magnetic field, which 
varies along the x direction. In order to study an smooth magnetic barrier we select a profile given by 

B = B Sech 2 (^) e z , (1) 

here e z is the unit vector normal to the graphene plane. This expression for the magnetic field presents several 
advantages: (i) We obtain an analytically solvable model, that allows us to analyze in detail the bound spectra, as 
well as the transmission through the magnetic barrier, (ii) In order to have conditions that are physically relevant to 
the study of graphene, we need a magnetic field B(x) that varies slowly on the scale of the graphene lattice spacing, 
a = 0.246 nm. Selecting a <§; d, we observe that both the half-width (A ps 3.25 c?) and the edge smearing length 
\(l/B)(dB/dx)\ ~ d of the magnetic barrier satisfy the required conditions. (Hi) As shown in Fig.((T]) the magnetic 
field in Eq. ([1]) provides a good approximation to the shape of the magnetic field barrier that is produced by a 
ferromagnetic film deposited on the top of a two dimensional system^!. Thus the present formalism could be useful 
to analyze a similar arrange of inhomogeneous magnetic barriers in graphene samples. 

The gauge can be selected in such a way that the vector potential is written as A = e y A(x) with 

A(x)=2B dtanh(^-A . (2) 

Notice that if we consider that the system is confined within a square box of area LxL, then in the d » L limit, 
the magnetic field can be consider homogeneous, and the vector potential reduces to the Landau gauge expression 
A = B o (0,x,0). 

On low energy scales the dynamics of quasiparticlcs in graphene is described by two independent 2+1 dimensional 
Dirac equations, the equations remain decoupled in the presence of smoothly varying magnetic field. The resulting 
time-independent Dirac equation describing low energy excitations around the K point in the Brillouin zone is written 
as 

Hip(x,y) = v F a ■ [p + eA(i)j ip(x,y) = Eip(x,y) , (3) 

here the Fermi velocity is vf ~ c/300, p = —ih\7 is the momentum operator and the isospin Pauli matrices Oi operate 
in the spinor tp(x, y) = (tpA , "0s) ; that represent the electron amplitude on two sites (A and B) in the unit cell of 
the graphene lattice. Taking into account the translational invariance along the y direction we seek solution of the 
form tpA = exp(ik y y)ip+ and ips = exp(ik y y)ip-. The Dirac equation yields the coupled equations 

A^+(x) - (-iJL-iW(x)\ i/>-(x), 

A0_(x) = (-i-^+iW{x)\i> + {x), (4) 



with the magnetic length defined as Ib = y/h/eBo, x = x/Ib, and A = Els/fiVF- The function W(x) is given by 

W(x) = l B k y + - h y ' = l B k y + — tanh (-) . (5) 



Combining the two equations in Q we obtain the decoupled equations 

dx 2 



H ± ^ ± (x)= l-_ + V±) i)±(x) = A 2 ip±(x) , 



where the effective potentials V± are given by 
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V±(x) = W 2 ± 



dW 
dx 



2d , / x 
— tanh — - 
l B \2d 



± sech 2 



(2d) 



(7) 



Both the shape and the depth of the effective potentials depend on the transverse momentum k y . As seen in Fig.©, 
depending on the values of k y and d, the effective potentials can assume the form of potential wells or steps. In next 
sections the bound state spectrum and scattering properties will be analyzed in detail. 

It is interesting to point out that the Dirac equation in the presence of an external magnetic fields posses a formal 
structure that can be analyzed within the formalism of supersymmetric quantum mechanics (SUSY-QM) 2 ^— . The 
potentials V+ and V- in are known as the super-partner potentials, they are obtained from the superpotential 
function W(x) by the relation in ([7]). The explicit expressions for V± for the gauge potential in are identified 
as the Rosen-Morse II potential o 28 i 29 and the corresponding Schrodingcr equations are exactly solvable. There are 
important property of SUSY-QM that relate the spectrum and cigenfunctions of the effective hamiltonians of H + and 
H- in Eq. (j6)). In particular, except from the ground state, H + and 7J_ have the same spectrum for A 2 . 

It is convenient to define the operators 



L± = -i—±iW(x) . 
dx 

In terms of these operators the relation between the upper and lower spinor components in @ simply read 



(8) 



Let us introduce the dimcnsionless variable 



tp-{x) = —L + tp + (x) . 



(9) 



1 + exp(x/d) ' 

that varies from to 1, as a; goes from 00 to —00. In the new variable equations ^ become 



(10) 



d 2 l-2£ d 
~dl 2+ W~0dt; 



A - [kyl B + 2d/l B {\ - 2£)] 2 ± 4£(1 - 



These equations have the asymptotic solutions: ip± 
where the asymptotic behavior is determined by 



£ 2 (l-£) 5 



0. 



^ for £ ^ (x 00) and ip± ~ (1 - £) a for ^ 1 (x 



(11) 

-00), 



(12) 



In order to obtain consistent solutions, we recall that besides solving the effective Schrodingcr equation in ([TTj) . the 
wave function components are interrelated by the Dirac equation via ©■ Then, one can consider the following options: 
(a) Equation (|11[) is solved for the lower component "0_ , the corresponding upper component ip+ is obtained from the 
first relation in ©. (b) Equation (jll[) is solved for the upper component ip + , and the lower component is obtained from 
the second relation in ©. We consider the first option, the second option gives equivalent solutions, except for the 
n = state. Taking into account the asymptotic behavior, we propose an ansatz of the form = £ p (l — 0^/(0) 

substituting in Eq. (|11[) we find that /(^) satisfies the Hypergeometric equation. Two linear independent solutions 
can be chosen as£i: /(£) = F (a, /3, 7; and /(0 = C 2p F (a - 7 + 1, /3 - 7 + 1, 2 - 7; £); where F (a, (3, 7; is 
the Hypergeometric function. However, the second solution has not the correct asymptotic behavior and has to be 
discarded. The corresponding upper component is obtained from the first equation in (|9]). The complete spinor 
solution is then given as 



Ad 



G(0 F (a, f3, 7; + 6(1 - Of F (a + 1, + 1, 7 + 1; 
F(a,P, r ,0 



(13) 
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where G(£) = Up — 2d 2 — £) — £(er — 2d 2 /l 2 B ) ~ k y d\ , C is the normalization constant and 



d\ 2 „ . ( d^ 2 



a = p + cr — 4: ( — J , /3 = p + CT + 4[ — ) +1, 7 = 2/7+1. (14) 



III. BOUND SATES 



Our aim is now to discuss the bound states spectrum. First we notice that the the electron-hole symmetry is 
preserved by the inhomogencous magnetic field. This follow from the fact that the Hamiltonian in (|3|) anticommutes 
with the (T3 matrix: {H, 03} = 0. Then, if ip is an eigenvector of H with eigenvalues E, a^ip is also an eigenvector 
with eigenvalue — E. 

The fact that the Hamiltonians H + and i/_ share the same eigenvalues, implies that the existence of bound states 
requires that both V + and V_ have the form of potential wells. This condition is obtained if the minima of V± are 
attained for a finite value of x, it is given by 

l 2 B \ky\<2d. (15) 

As seen in Fig.([5Ji), when the previous condition holds both V+ and V- have the form of asymmetric potential wells 
around the guiding center position 



x c = 2darctanh ( ^ j . (16) 

Bounds states are found if the wells are sufficiently deep. In the limit in which the transverse momentum vanishes 
(k y = 0), the effective potentials reduce to symmetric wells centered at the maximum of B{x). For values of k y in 
which the condition (|15|) is not obeyed at least one of V± take the form of a potential step Fig.([5]3), and bounds states 
are not supported. 

The solution in Eq. (|13p leads to a divergence in £ — > l(x — > —00) except for a or j3 being a negative integer. 
Letting a = —n, and utilizing Eq. (|T4"]) we obtain the energy spectrum, that can be conveniently written as 



E n .k 



± 



2hev 2 F B Q 



\ 



(m B y 

8d 2 



1 - 



3 ky/2d 



1 - n{l B /2df 



(17) 



The index n take the values n = 0, 1, 2, 3. 



Both the allowed values of n and and k„ are restricted in order to 



satisfy the square integrability condition, explicitly they yield 



/ d\ 2 d Ud 2 -nl 2 } 2 , N 

The first condition determines the highest bound state supported for a given width of the barrier and eliminates the 
possible singularity in (|17p . Whereas, the second condition determines the allowed values of k y for a given n and is 
related to the fact that the electron group velocity is limited by the free velocity in graphene, i.e. \dE ni k y /dp y \ < vf- 
The two conditions guarantee that E n ^ < (fi-VF/^s) rnin {V±(±oo)}, so the electron does not escape towards ±00; 
cquivalently the coefficients in (|12[) that determine the wave function asymptotic behavior satisfy p, a > 0. 

The wave function for the zero energy level takes a simple form that can be obtained by solving the first equation 
in ((4]) when A = 0, it reads 



i' = Ce lk * y e- k » x \Sech(^\ 



i(d/h 



,2/0 



1 



(19) 



Whereas, for other values of n the Hypcrgcomctric functions become Jacobi polynomials with energy dependent 
indices, the wave function is given by 
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fa M{z) p£- x <- n+ ^\z) - (1 - z 2 ) f£ P^- 1 ' 1 -"^-^^) 
tjj = Ce ikyV (l + z) p (1 — z) a | |, (20) 

p h,-i,-„+0- 7 )^ 

where 2 = 2£ - 1 = tanh(x/2) and M(z) = (p - 2d 2 /l 2 B )(l - z) - [a - 2d 2 /l%)(l + z) - 2k y d. 

The dispersion relation in (|17|) shows that the inhomogeneity of J5 lifts the degeneracy for every quantum level n 
and gives rise to a fc y -dependent dispersion relation, which leads to a drift velocity along the y axis. This is valid, 
except for the n = level that has zero energy, independent on the magnetic field for all values of k y . The energy 
spectrum for d = 1.51b is displayed in Fig. ([3]) as a function of k y . For the selected values we have n max = 9, and each 
level n results in a band as the transverse momentum sweeps from k y — to the maximum value k y max in l|18p . 

The previous results apply when the width of the barrier is small in comparison with the system dimensions. We 
now analyze the behavior of the spectrum as we change from a narrow to a broad barrier, considering that the system 
is confined within a square box of length L x L. For a narrow barrier the values of k y are limited by the second 
equation in (|18[) . Instead, when the size of the barrier is comparable to L, the number of allowed states is limited by 
those that can be accommodated in the square box. Assuming periodic conditions for the wave function in (|20l) along 
the y— direction, yield k y = 2ir j / L with j an integer. But according to Eq. (1161) k y also determines the center position 
x c of the electron, hence x c < L/2 and the number of quantum states is given by N = j max = (2dL/irl B ) tanh (L/Ad), 
whereas the momentum limit imposed by the size of the system reads k y max = 2(d/l B ) tanh (L/Ad). It is interesting 
to notice, that similarly to the homogeneous case, the degeneracy can be written as TV = $/$o, where the magnetic 
flux produced by the field in ([1]) through the sample is $ = AB^dLtavAi (L/Ad) and <t>o = h/e is the elementary fluxon. 
In order to track the degeneracy evolution as the barrier is modified from narrow to broad as compared to L, we 
define 



Hy,max ^ yi max ""y,max 



as a cut for the transverse momentum. Notice that K Vt ma x interpolates between k ymax valid for a narrow barrier, and 
ky max valid when d > L. Fig. ((4]) shows the resulting energy spectrum, the dark zones are the allowed energy values. 
For every level n the transverse momentum varies between k y = and k y = K Vt max . The restriction on the level index 
n given by the first equation in (|18[) . translates into the following equation for the separatrix E = 2hvpd/l 2 B (dashed 
line), energies to the left of this line are not allowed. For small value of d/ls a few bands and gaps can be identified 
for the first values of n, as the energy is increased we observe a continuous energy region up to the maximum allowed 
value n max . As the size of the barrier increases a larger number of energy gaps appear. When the size of the magnetic 
barrier becomes comparable to L the width of the bands decrease. Finally in the homogeneous B limit, d^ L, the 
energy eigenvalues reduce to the relativistic massless Landau levels observed in graphene: E n ^ y = i\/2hevpBo n, 
independently of k y , whereas the degeneracy reduces to the well known result N = $/3>0: with $ = B L 2 . 



IV. SCATTERING 



We now consider the scattering regime. A plane wave incident from x — > — oo , (£ — > 1) propagates at an angle <$> 
with respect to the x axis. Taking into account the gauge selection in ([2]) the incoming momenta are parametrized as 



fcx = A cos , Ib k y = As'mcj) + — . (22) 

Ib 

The transmitted wave has a longitudinal momentum k' x = C os <fi' where <p' is the refracted angle. The conservation 
of k y gives the relation between the incident and refracted angle as follows 



sin (j)' — sin <f) + - — — , (23) 
Ib A 

whereas energy conservation allow us to relate the transmitted and incident longitudinal momenta as 
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l B k' x = \J(lBk x ) 2 - \J^J -8k y d. (24) 

Eq. (|23p implies that for a critical angle <f> c = arcsin(l — 4d/l B A) no transmission is possible. Furthermore, when 
the following condition applies 

El ft d 

T-^<2 r , (25) 

the transmission vanishes regardless of the incident angle <f>. The condition in (|2"5"|) establish that states with an 
average cyclotron radius (in the barrier region) smaller that 2d will bend by the magnetic field and are completely 
reflected, a similar condition was obtained in the case of a square-well magnetic barrier—. Comparing the equations 
(f2"2"|) with (|12p . we observe that the longitudinal momenta k x and k' x are related to the asymptotic coefficients p and 
a as follows: a = —ik x d and p = d \/—k' x . It is verified that if condition (|25|) holds p is real, instead when (|25|) is not 
valid we replace p = —ik' x d. Utilizing the properties of the Hypergeometric functions it is verified that in the limit 
x —> oo (£ — > 0) the wave function in (|T^|) yields the correct asymptotic expression 



4rf 



2 



where s = sgnE. The asymptotic value of the wave function for x — > — oo (£ — > 1), is obtained using the linear 
transformation formulas^ that relate F(a, ft, 7, £) with Hypergeometric functions evaluated at 1 — £, to obtain 



From this equation the reflection coefficient i? is obtained as 



i? = 



r ( P - it x d + (2d/i B ) 2 + 1) r ( P - it x d - (2d/i B ) 2 ) 



T(p + ik x d + (2d/l B ) 2 + 1) r (p + ik x d - (2d/l B ) 2 ) 



(28) 



Under condition ([25]) p is real, thus R = 1 and, as expected, the transmission coefficient vanishes. Instead, if the 
condition (|25|) is not obeyed, p is substituted by p = —idk' x , and the transmission probability T = 1 — R can be 
written in a simple closed form as 



T 



sinh [2irl B k x ] sinh [27r/sfc^.] 



sin 



4tt 



(29) 



sinh 2 [ivl B {k x + k' x 



A plot for the transmission coefficient as a function of the incidence angle is shown in Fig.© for a fixed energy 
and several values of d. The qualitative behavior in this plot is similar to those obtained for the case of a square 
well magnetic barrier—. However an advantage of the the result in (f29| is that it allow us to identify the resonant 
conditions at which the barrier becomes transparent (T ~ I), it is given by 



d=} y l B .7 = 1,2,3, (30) 

In this case the barrier acts as an asymmetric filter, it behaves as perfectly transparent for angles in the region 
— 7r/2 < cj) < cj) c Fig.([l)k.). In particular for energy values slightly above the threshold condition in (|2"5"|) . A = d/l B + e, 
with e <C 1, the width of the transparency region can be very narrow. 
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On the other hand, under the condition 



d = 




3 = 0,1,2,3 



(31) 



2 



the value of T is reduced, in particular for values of the energy slightly above the threshold condition in (|25j). 
A = d/ls + e, with e <C 1, the transmission coefficient is strongly reduced, Fig. (J&h) . Another way to contrast these 
results is shown in Fig.© where a contour plot for the transmission coefficient T is shown. A can be seen in Fig.([7ta) 
corresponding to the resonant condition in Eq. ([30]) there is a wider region for the transmission that is not present in 
Fig-Gib) corresponding to the minima condition in Eq. (|31|) . 



In conclusion, we consider the dynamics of carriers in graphene subjected to an inhomogeneous hyperbolic magnetic 
field. The corresponding Dirac equations was analyzed within the formalism of supersymmetric quantum mechanics. 
We found compact analytical solutions for the energy eigenvalues and eigenfunctions for electrons and holes. The 
dispersion relation in (|17|) shows that the inhomogeneity of B lifts the degeneracy for every quantum level n and gives 
rise to a fcy-dependent dispersion relation, which leads to a drift velocity along the y axis. This is valid, except for 
the n = level that has zero energy, independent on the magnetic field for all values of k y . For a narrow barrier 
the spectra displays a series of bands separated by gaps, as the width of the barrier increases we can track the levels 
evolution into the degenerated Landau levels. In the scattering regime a simple analytical formula is obtained for 
the transmission coefficient, this result allow us to identify the resonant conditions at which the barrier becomes 
transparent. We expect that the results obtained in this work will be useful in order to analyze the confinement by 
magnetic barriers in graphene samples. In a future work we plan to address the problem of calculating the longitudinal 
conductivity as well as the Hall conductivity for the case of graphene under inhomogeneous magnetic fields. 
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FIG. 1: (color online) Magnetic field profile for the hyperbolic expression in Eq. fTJ (solid line) as compared to the profile 
B z (x) — — (/UqM/47t) ln[(s 2 + k 2 )/(x 2 + (I + /i) 2 )],— produced by a ferromagnetic strip with a magnetization fXoM, thickness 
I and separated a distance h from the two dimensional system (dashed line). The parameters are selected in such a way that 
the peak and half width of both profiles have the same value. 




FIG. 2: (color online) Plots of the effective potentials V- (solid lines) and V+ (dashed lines) versus x/Ib- (a) The parameters 
kyls = 1, d/ls = 2 satisfy the condition in (|15p . The A 2 (k y ) eigenvalues are represented by the dotted lines, (b) kyls = 1, 
d/ls = 0.4, notice that in this case the condition in (|15|l is not satisfied, and therefore bounds states are not supported. 




FIG. 3: (color online) Bound states energy spectrum E n (k y ) as a function of k y for a magnetic barrier with d — 1.5 Ib- According 
to (|18|l Umax = 9. For every n— level the allowed values of k y are delimited by the free-electron spectrum E = ±hvpk y (dashed 
lines) . 
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FIG. 4: (color online) Energy spectrum as a function of d for a magnetic barrier when the system is confined by a square box 
of area L x L with L = 50/s. The dark zones are the allowed energy values. The dashed line EIb/^vf = 2d/ls delimitates the 
minimum barrier width that supports bound states. For every level n the transverse momentum varies between k y — and 




FIG. 5: (color online) Angular dependence of the transmission coefficient T through a barrier with energy A = EIb/Uvf = 3 
and different values for d: d — \AIb (continuos line), d — Ib (dashed line), and d — 0.51b (dotted line). 
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FIG. 6: (color online) (a) Angular dependence of the transmission coefficient T through a barrier with d = Zb/2 (resonant 
condition in Eq. (|30[) with j = 1), and the following values of the energy: A = 1.01 continuos line, A = 1.03 dashed line, and 
A = 1.1 dotted line, (b) As in (a) for d = Zs/Vo (Eq. (|31[) with j = 0) and the following values of the energy: A = 0.75 
(continuos line), A = 0.8 (dashed line), and A = 0.85 (dotted line). 




FIG. 7: (color online) a) Contour plot of the transmission coefficient T through a barrier with d = Ib/2 (resonant condition in 
Eq. (O with j = 1), (b) As in (a) for d = l B /V8 (Eq. J3T) with j = 0). 



